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The equilibrium structures and electronic excitation spectra of the Ir(III) photosensitizer 
[Ir(ppy) 2 (bpy)]^ bound to medium-sized silver clusters Ag^ (n =19, 20) are investigated using 
time-dependent density functional theory. The long-range corrected LC-BLYP approach is used 
with a system-specific range separation parameter. The weak physisorption of the hybrid complexes 
yield only small changes in the broadened absorption spectra of the hybrid system as compared 
with its constituents. However, the density of states as well as the fine structure of the spectra is 
strongly modified upon complexation. It is shown that the standard range separation parameter 
(0.47 bohr“^) cannot predict these properties correctly and the optimized value of 0.16 bohr“^ 
should be used instead. 


I. INTRODUCTION 


The combination of metal nanoparticles or nanowires 
with various organic adsorbates such as dyes, peptides, 
and J-aggregates is an active area of research with ap¬ 
plications in bio-sensoring, catalysis, and medicine 
[5]. Numerous theoretical and experimental studies of 
nanoparticle-organic hybrid systems have been reported 
(for reviews, see, e.g.. Refs. [6]-[9]). These hybrid systems 
exhibits new composition dependent properties, which 
differ from those of the separate constituents This in¬ 
cludes the enhancement or quenching of fluorescence, ab¬ 
sorption, and Raman scattering due to surface and plas- 
mon resonances, the broadening of the absorption range 
to yield an antenna effect, and modified “redox” proper¬ 
ties; see, e.g.. Refs. [IMl]. 

In the present communication we discuss results on 
the electronic properties of a model system comprised 
of the [Ir(ppy) 2 (bpy)]^ photosensitizer (IrPS) shown in 
Fig.[^ bound to medium-sized silver clusters Ag^ (n=19, 
20). Our choice of the system is motivated by the use 
of heteroleptic Ir(III) complexes as photosensitizers in a 
photocatalytic system for water splitting introduced by 
Beller et al. m The original homogeneous system in¬ 
cludes triethylamine as a sacrificial reductant and a series 
of iron carbonyls as water reduction catalysts. The ab¬ 
sorption spectrum of the IrPS overlaps with the sun’s 
spectrum only in its long wave length tail between 300 
and 450 nm m- Therefore, extending the range of ab¬ 
sorption further into the visible by coupling the IrPS to 
silver nanoparticles might provide a means to enhance 
the overall efficiency of this photocatalytic system. 

Because of size of the systems under study, density 
functional theory (DPT) [17] and its extension into the 
time-dependent (TD) domain in the linear response limit 
m is the natural choice for studying ground and ex- 
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Figure 1: Chemical formula of the photosensitizer 

[Ir(ppy) 2 (bpy)]+ (IrPS). 


cited state properties. However, the correct description 
of long-range charge-transfer (CT) properties, which is 
of mandatory for hybrid systems, can not be achieved 
applying conventional DFT functionals [TM23] . By in¬ 
troducing exact Hartree-Fock exchange in long-ranges 
corrected schemes as LC-BLYP [241126] or CAM-B3LYP 
m, the correct asymptotic behavior and more balanced 
description for CT states is obtained. In our previous 
study on bare IrPS m, the LC-BLYP approach was 
shown to give excitation energies in good agreement with 
experimental data and results of CASSCF/CASPT2 cal¬ 
culations. 

The binding of IrPS to small model silver clusters (one 
to six silver atoms) has been studied using the LC-BLYP 
approach in Ref. [28]. Already for this size of the hybrid 
system pronounced changes in the electronic absorption 
spectrum had been observed. Moreover, the electronic 
properties were found to be rather sensitive to the num¬ 
ber of silver atoms. Binding energies and localization of 
HOMO and LUMO orbitals are found to oscillate with 
the number of silver atoms in an even/odd like fashion. 
Whether such a behavior propagates to larger clusters 
and how it influences the photophysical properties is the 
topic of the present study. 

The paper is organized as follows. In Sec.[n]we give the 
computational details putting emphasis on the tuning of 
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the range-separation functional in LC-BLYP. The results 
are presented in Sec. m including electronic ground state 
structures and electronic absorption spectra. A summary 
is provided in Sec. EYl 


II. COMPUTATIONAL DETAILS 

In the long-range separation approach LC-BLYP a pa¬ 
rameter uj is introduced, which defines the separation of 
the Coulomb operator into long-range and short-range 
parts, with the long-range part then being described by 
the exact exchange integral. In a number of publica¬ 
tions, it was demonstrated that an appropriate tuning of 
UJ leads to a significant improvement of fundamental and 
optical gaps, CT and Rydberg excitation energies as well 
as ionization potentials (IPs) [29H32]. A systematic pro¬ 
cedure for the determination of an optimal uj has been 
suggested in Refs. [33II35] . It is based on finding that cj, 
which minimizes the following functional 

J{ij) = J]sf{ij) + 

= kSoMoW+IP^WI 

+ noMoiN + l) + lP^{N + l)\. ( 1 ) 

In other words, the optimized uj will provide a compro¬ 
mise for fulfilling Koopmans’ theorem simultaneously for 
systems with N (cation complex) and Y -|- 1 (neutral 
complex) electrons. In Eq. Q, ^homo 
HOMO energy and the IP, respectively. 

In Ref. [36] we have studied the present systems 
and obtained optimized range-separation parameters for 
the bare IrPS, silver clusters Ag^, and hybrid systems 
IrPS—Ag^ (n = 2, 10, 20). Here, we give a more detailed 
account on the case of IrPS—Ag^^Q. In short, the IPs and 
HOMO energies of neutral species, anions, and cations 
have been obtained from single-point energy calculations 
at the optimized geometry of the electronic ground state. 
Similar to the small systems [28|, the geometry opti¬ 
mization was carried out without symmetry constraints. 
Initial geometries of IrPS for further optimization were 
taken from Ref. [16] (assuming C 2 point symmetry |37|). 
The initial structures of Ag^^Q were taken from Ref. [38], 
those of Ag^g and Ag 2 o from Ref. [39], except tetrahe¬ 
dral Ag 2 o which had been studied before in Refs. [40ti42] . 
Binding energies have been obtained including the coun¬ 
terpoise method to correct for the basis set superposition 
error (BSSE) [43l[44] . 

Vertical excitation spectra have been calculated by 
means of the TDDET approach. The number of included 
transitions has been 450 and 650 for XXa and IXXo, re¬ 
spectively. A Lorentzian broadening (0.4 eV) has been 
added to the stick-spectra. Spin-forbidden transitions 
were not included in the present TDDET investigation 
although they could have notable intensity due to the 
high spin-orbit coupling constant of Ir. All calculations 
were performed with the GaussianOO program package 



CO, bohP 

Figure 2: Optimization of range-separation parameter for 
the case of IrPS—Ag^g. The functions Jn{^) and 
are defined in Eq. 0 and represent the conditions for fulfill¬ 
ing Koopmans’ theorem for the cationic and neutral species, 
respectively. 


m using the LANL2DZ effective core potential basis set 
for Ir and Ag and the 6-31G(d) basis for all other atoms. 


III. RESULTS AND DISCUSSION 
A. Range-separation Parameter Optimization 

In Eig. the procedure for tuning uj is demonstrated 
for the case of IrPS—Ag^^o- ^ minimum of J{uj) at 0.16 
bohr“^ was found. Around 0.63 bohr“^ an additional 
higher energetic minimum occurs. This can be traced to 
the change of the character of the lowest adiabatic state 
of the cationic species, where the order of HOMO and 
HOMO-1 orbitals are interchanged. Here, the unpaired 
electron initially localized on the 7r*(bpy) orbital of the 
reduced doublet IrPS is transferred to the cr(Ag;Lo) 

bital [in]- 

The obtained optimal uj value is substantially smaller 
than values of 0.33bohr~^ and 0.47bohr“^ (shown as 
vertical lines in Eig. implemented in common quan¬ 
tum chemical packages. This is not surprising since these 
(standard) values were determined for test sets of di¬ 
atomics and small molecules [251147]. 

Since the inverse of uj reflects a characteristic distance 
for switching between short- and long-range parts of 
the exchange contribution, optimal uj values were shown 
to decrease with increasing system size and conjugation 
length EniisniisaissiisHsn]- This size-dependence, how¬ 
ever, is not monotonous and there is a strong dependence 
on the electronic structure [32]. Optimal values of uj near 
0.2 bohr~^ have been reported for various systems with 
characteristic electronic radii comparable to the present 
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Figure 3: Optimized structures of IrPS—Ag^^g and 

IrPS—Ag 2 o studied in this paper. The conformers labeled 
IXXe and IXXo are not distinguishable at the scale of the 
image. 

ones, e.g., 0.214 bohr“^ was found for pentacene-C 0 o |5H : 
for other examples see Refs. [saissisH]. 

As shown in Ref. [36] there is only a slight variation 
of uj when going from IrPS, via Ag^ (n = 10 , 20 ) to 
IrPS—Ag]^Q 20 - Thus, in the following we provide results 
on excited states obtained using uj = 0.16bohr“^ for all 
compounds. It turned out the ground state geometries 
of the present systems are little affect by the choice of uj. 
The corresponding tests have been performed for IrPS, 
Ag 2 , IrPS—Ag 2 , and IrPS—Ag^g. Therefore, we have 
used the the standard co ’2 = 0.47bohr“^ for the ground 
state optimization. Selected TDDFT calculations are 
also performed with the latter value in order to study 
the influence of uj on the absorption spectrum. 


B. Ground electronic state 

In our previous study of IrPS bound to small (n < 6 ) 
silver clusters [28], we found that configurations in which 
Ag^ is situated in the cavities between ligands are the 
lowest in energy. Note that in all cases the interactions 
are “weak” and no chemical bonds are formed. Here, 
we extended our investigation to IrPS—Ag^ geometries 
with 19 and 20 silver atoms, however, focussing mostly 
on those structures where the cluster is located in the 
“ppy-ppy” cavity as shown in Fig. 

The differences in geometry of IrPS upon complexation 
with a silver cluster are only minor (maximal changes are 
0.09 A and 4.6°); for an overview of geometries see Table 

The binding energies of the studied systems are in the 
range 1.5 — 2.5 kJ/mol, which is about 10 times smaller 
than that of silver atoms in large nanoparticles [52] and 
smaller than in analogous systems with n = 1 — 6 . This 
confirms the trend of a decrease of binding energies with 
the increase of cluster size found previously [28]. Such 


a behavior can be explained by the fact that no stable 
complexes are formed, i.e. the interaction belongs to ph- 
ysisorption case and might be considered as a weak at¬ 
traction between the ligands of IrPS and nearest atoms of 
silver cluster. The total binding energies (not divided by 
number of silver atoms) are 30-100 kJ/mol for all systems 
disregarding the cluster size. Accounting for the BSSE 
decreases the binding energy by factor of about two for 
all cases. The change of uj from the standard to the opti¬ 
mized value has no uniform impact on binding energies. 
Changes of a factor of about two in both directions were 
observed. Note that the order of stability is sensitive to 
the value of range-separation parameter used. For clar¬ 
ity, we provide here only results of calculations with the 
optimal UJ value. For the series of XX stuctures, the XXa 
conformer is the lowest in energy, while XXa’ and XXb 
are by 30.0 and 4.1 kJ/mol higher, respectively. For the 
IXX species, IXXe is more stable (AE = 38.2kJ/mol). 
Nevertheless, in the following we also discuss the IXXo 
structure because of its unique electronic structure in the 
ground electronic state. 

The representative frontier orbitals (LC-BLYP, uj = 
0.16 bohr“^) are given in Eig. for the cases of IXXo 
and XXa. If the standard value of uj is applied then the 
HOMO-LUMO gap increases by about 1 eV for all cases. 
The shift is mainly due to the LUMO orbitals while the 
nature and order of orbitals remains the same. 

The orbital structure differs for systems containing odd 
and even number of silver atoms. Eor an even number, 
the HOMO orbital is localized on the silver fragment 
(^Ag) LUMO is of 7 rjf(bpy) type. Eor an odd 

number of silver atoms this situation is reversed: the un¬ 
paired electron is situated on the 7rjf(bpy) orbital as in 
the reduced IrPS, while the LUMO is localized on the 
silver fragment. A similar situation with one exception 
for IrPS—Ag 5 was found for smaller clusters [28] as well 
as for the much larger cluster IrPS—Ag 92 |53| . 

Eor 19 atoms two cases were found which slightly differ 
in geometry: IXXo (“o” means “odd”) follows the orbital 
rules just given and IXXe (“e” stands for “even”) is an 
exception where the unpaired electron is localized on the 
silver fragment. The “odd” and “even” groups can be also 
seen from the comparison of geometries, see Table [Tj the 
distances Ir-N^p^ and r(C-C) 0 py are shorter for the “odd” 
group and the IrPS in its reduced form (IrPS^). 


C. Absorption spectra 

The absorption spectra of IXXo and XXa are presented 
in Eig. The standard (0.47 bohr“^) and optimized uj 
values were applied for comparison. Because of the num¬ 
ber of states considered and of the complex character of 
most states, we distinguish only between two types of 
orbitals: those localized on silver or on dye fragments, 
but do not consider the particular character of orbitals, 
e.g., 7 r(ppy), 7 r(bpy), or (i(Ir) for IrPS. Eor each state 
different contributions have been summed up and the 
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Table I: Selected equilibrium geometric parameters of all IrPS- • • Ag^ compared to IrPS [16] and IrPS° [46] in the 
ground electronic states. For notation of structures see Fig. Bond lengths in A, angles in degrees. 


Parameter ^ 

IrPS 

IrPS"" 

IXXo 

IXXe 

XXa 

XXa' 

XXb 

r(Ir—Cppy) 

2.007 

2.010 

2.019 

2.005 

2.007 

2.001 

2.002 

r(Ir-Nppy) 

2.055 

2.051 

2.051 

2.059 

2.058 

2.056 

2.057 

r(Ir-Nbpy) 

2.163 

2.134 

2.116 

2.158 

2.158 

2.160 

2.162 

r(C-C)ppy 

1.463 

2.464 

1.470 

1.463 

1.466 

1.457 

1.461 

r(C-C)bpy 

1.485 

2.409 

1.408 

1.482 

1.485 

1.481 

1.482 

r(Ir-Ag) ^ 

- 

- 

3.563 

3.919 

4.042 

5.101 

5.108 

ZCIrC 

89.3 

89.1 

93.9 

92.9 

90.7 

89.8 

90.4 

ZNppyIrNppy 

172.1 

174.5 

174.8 

171.8 

171.8 

172.8 

172.2 

dCIrCN 

92.4 

96.1 

95.8 

93.2 

93.3 

95.0 

94.6 


^ For parameters including ppy, only those ppy fragment that is closer to the silver cluster is regarded. The pa¬ 
rameters of second ppy fragment differ by no more than 0.04 ; ^ Distance between Ir and closest silver atom. 


LUMO _ LUMO 


a*(Ag) - ^ TT,*(bpy) 



- £orb. eV 

IXXo XXa 


Figure 4: Molecular orbital energies and frontier orbitals for 
IXXo and XXa systems. 

corresponding dominating character has been assigned. 
Most of transitions have very complex assignment and 
are marked as “mixed”. In Fig. |^the spectra are shown 
as sticks and with a Lorentzian broadening. 

All spectra have a clear maximum composed mainly by 
silver-localized transitions, which, for pure cluster spec¬ 
tra, would be called “plasmonic” excitations (collective 
excitation of electrons). Similar to the cases of small 
silver clusters [28], the hybrid systems support a new 
type of long-range intermolecular CT, corresponding to 
IrPS ^ Ag^ or Ag^^ -> IrPS. These intermolecular CT 
transitions (denoted as IM to emphasize the their long- 
range nature) are marked with red bars in Fig. The 
IM transitions have very different intensities, from almost 
zero to that comparable with silver plasmon-like transi¬ 
tions. Many transitions can not be attributed to pure IM 
type, but have a notable contribution of intramolecular 
CT character. These mixed transitions are marked with 
black bars in Fig.[^ In the lower panel of Fig.[^ electron 
density difference plots for transitions of IM and cluster 
localized types are given. 

Comparing the spectra of XXa and IXXo with different 


uj values, three important points should be highlighted. 
First, their maxima are shifted to the red by 0.4 eV for 
XXa and by 0.1 eV for IXXo. The positions of particular 
bands can hardly be compared as the assignment is very 
complex. Second, the density of states for the optimized 
UJ is almost twice higher than for the standard uo. Finally, 
the distribution of oscillator strength depends on uj as 
well. This is particularly visible for the XXa system, 
where the optimized uo yields a double peak spectrum, 
which is in contrast to IXXo. For the standard uo XXa 
and IXXo have essentially the same line shape. 

The changes of absorption spectra upon complexation 
are of general importance for possible photochemical and 
photophysical applications. Because of complicated as¬ 
signment of single bands, the direct comparison of partic¬ 
ular transitions is hindered and only the overall shapes 
can be compared. In Fig. the broadened spectra for 
IXXo and XXa are compared to those of the pure con¬ 
stituents. For IXXo, we provide the spectra of reduced 
IrPS^ and Ag^g according to the ground electronic state 
structure, whereas for XXa the results for IrPS and Ag 2 o 
are given. 

The broadening parameter was chosen in such a way 
as to reproduce the approximate width of experimental 
plasmonic band for silver clusters, which is about 50 nm 
[5l [55]. However, different types of transitions (local 
and CT), even in case of pure IrPS, could have a dif¬ 
ferent broadening, which is not considered due to the 
lack of detailed information. For the case of smeared 
fine structures of spectra, we expect our results to be 
only slightly dependent on the variation of line shape 
and small changes in line widths of transitions of partic¬ 
ular types. When discussing the changes in the broad¬ 
ened spectra in the following, one should keep in mind 
that the individual transitions and thus the unresolved 
fine structure of the spectrum is strongly modified upon 
complexation. 

If we compare the absorption spectra of the hybrid 
systems with those of their constituents the following 
observations can be made: The maximum of the plas¬ 
monic band in the spectrum of IXXo is shifted by 0.2 
eV to the blue and its width increases slightly. For XXa 
the maximum remains almost at the same position upon 
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Figure 5: Upper panel: Selected vertical electronic singlet-singlet (for XXa) and doublet-doublet (for IXXo) TDDFT spectra 
calculated with LC-BLYP with standard and optimized u) parameter. Associated spectra are broadened with Lorentzians of 
width 0.4 eV. Lower panel: Electron density differences (contour value 0.0008) for selected transitions, which are marked with 
asterisks in upper panel. Red and blue colours correspond to particle and hole densities, respectively. 


complexation with IrPS. The position and shape of this 
broad band is very similar for unbound and bound clus¬ 
ters. Some additional features in the spectra of XXa in 
the region from 4.5 to 5 eV could be attributed to the in¬ 
fluence of IrPS through the manifold of new states with 
non-negligible IM CT character. 

Figurej^also contains the sum of spectra of the isolated 
dye and the cluster. It almost does not differ from the 
spectrum of the interacting hybrid system for XXa, and 
it is shifted only by about 0.2 eV to the blue for IXXo. 
But the density of states increases and most transitions 
have at least an admixture of IM character. Therefore, 
the spectra of hybrid systems cannot be considered as 


perturbed spectra of metal clusters. Consequently, the 
photodynamics of the hybrid systems might strongly dif¬ 
fer from those of pure IrPS. 


IV. CONCLUSIONS 

The interaction of an Ir(III)-based photosensitizer with 
silver clusters containing 19 or 20 atoms has been stud¬ 
ied by means of TDDFT with a properly chosen range- 
separation parameter in the LC-BLYP functional. Al¬ 
though the interaction between dye and metal particle 
belongs to the case of weak physisorption, the properties 
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tion spectra of pure IrPS, Ag^, their sum, and combined sys¬ 
tems IrPS—Ag^ (n = 19, 20). The width of the Lorentzian 
broadening has been 0.4 eV. Note that for IXXo the reduced 
IrPS and the oxidised Agig have been used as reference, in 
accord with the ground state orbital structure of this system. 

of the hybrid systems differ notably from those of its con¬ 
stituents. This concerns the appearance of long-range in- 
termolecnlar CT states, which might lead to a completely 
different photophysical and photochemical behaviour, in¬ 
cluding long-lived charge separation favourable for fur¬ 
ther reactions. The low resolution electronic excitation 
spectra are close to the sum of the spectra of the non¬ 
interacting constituents. However, the underlying den¬ 


sity of electronic states and thus the fine structure of the 
spectra is strongly modified in the hybrid system, with 
the details depending on the number of silver atoms. The 
changes in the adsorption spectrum upon forming the hy¬ 
brid system are mostly observed the region from 3 to 4.5 
eV, but they are masked by the broad plasmonic type 
absorption band of the cluster. 

If we extrapolate the obtained theoretical results to 
experimental observations the following considerations 
should be taken into account. First, metal clusters used 
for spectroscopic and catalytic studies are produced ei¬ 
ther chemically by reduction in solutions [56l l57] or with 
laser techniques, e.g., cluster beam generation utilizing 
arc discharge, magnetron sputtering, or laser vaporiza¬ 
tion (for review see, e.g.. Ref. [58]). The resulting clusters 
are substantially larger compared to those used in our 
theoretical model and have a certain size distributions. 
Second, the clusters are produced either in solution or de¬ 
posited on a plate. The deposited nanoclusters have the 
advantage of being protected from fast aggregation and 
thus are often applied in spectroscopic studies and also 
in catalysis providing more stable and long-living cata¬ 
lysts. The experimental maximum of plasmonic band for 
silver particles is in the range 375-430 nm for an average 
diameter from 2 to 50 nm [S1HS71ISS] and produced in 
aqueous solution or embedded in solid substrate (Si 02 , 
Cr 203 , and MgF 2 ). The maximum of this peak shifts 
to the red with increasing the average cluster size and 
with increasing the dielectric constant of the surround¬ 
ing medium. Consequently, the spectra of IrPS bound to 
silver clusters relevant in those experiments should be av¬ 
eraged over different cluster sizes and red-shifted as com¬ 
pared to those shown in Fig. It can be expected that 
the differences between odd and even number of atoms 
in cluster will not be seen in this case. Further, the in¬ 
teraction between the silver clusters and the IrPS will be 
weakened due to the increase of the spectral separation. 
The consequences for the properties of the hybrid system 
still need to be explored. 
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